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ABSTRACT 

We report the earliest detection of an extremely bright optical afterglow of 
the gamma-ray burst (GRB) 030329 using a 30cm-telescope at Tokyo Institute 
of Technology (Tokyo, JAPAN). Our observation started 67 minutes after the 
burst, and continued for succeeding two nights until the afterglow faded below 
the sensitivity limit of the telescope (approximately 18 mag). Combining our 
data with those reported in GCN Circulars, we find that the early afterglow light 
curve of the first half day is described by a broken power-law (oc t~ a ) function 
with indices a x = 0.88 ± 0.01 (0.047< t < t hl days), a 2 = 1.18 ± 0.01 (t bl < t < 
tb2 days), and «3 = 1.81±0.04 (^2 < t < 1.2 days), where tbi ~ 0.26 days and ib2 
~ 0.54 days, respectively. The change of the power-law index at the first break 
at t ~ 0.26 days is consistent with that expected from a "cooling-break" when 
the cooling frequency crossed the optical band. If the interpretation is correct, 
the decay index before the cooling-break implies a uniform ISM environment. 

Subject headings: gamma-rays: bursts - radiation mechanisms: non-thermal 



1. Introduction 



The overall behavior of GRB afterglows has been successfully explained by the standard 
fireball model (e.g. Piran et al. 1999). It is expected that the deviation from a simple power- 
law in the afterglow light curve will provide wealth of information on the environment and 
physical parameters of GRBs. For example, observations of early GRB afterglows confirmed 
that "breaks" exist in the light curves of a number of GRBs. Such breaks may be understood 
in the framework of the standard fireball model either as a "jet-break" where the bulk 
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Lorentz factor of the relativistic jet decrease to the inverse of jet opening angle (Sari, Piran 
& Halpern 1999, Rhoads 1999), or the "cooling-break" where the high energy electrons start 
to lose most of their energy rapidly by synchrotron emission which emit the observed optical 
photons (Sari, Piran & Nakar 1998). While there are number of convincing cases for jet- 
breaks, identification of a cooling-break in the optical light curve has been difficult since 
it requires detection of a subtle change in the power-law index (Act ~ 0.25). In order to 
study afterglow light curves in detail, continuous coverage and a high signal-to-noise ratio 
are required. 

The situation has been dramatically improved since the advent of High Energy Transient 
Explorer-2 (HETE-2). HETE-2 can determine the positions of GRBs onboard, and notify 
ground-based observers of the GRB coordinates within 1 — 10 minutes after the burst (Ricker 
et al. 2002). For example, the locations of GRB 021004 (e.g., Fox 2002) and GRB 021211 
(e.g., Fox & Price 2002) were disseminated within less than a minute after the bursts, which 
prompted detailed studies of GRB afterglows in the very initial phases while they are bright. 

GRB 030329 was detected by the HETE-2 satellite on 29 March 2003 at 11:37:14.7 UT. 
Position was determined by the ground analysis, and the location was reported to GCN at 
73 minutes after the burst (Vanderspek et al. 2003). A very bright (R ~ 13 mag) optical 
transient (OT) was reported at a = 10 /l 44 m 50 s .O, 5 = +21°31'17.8"( J2000.0) (Peterson & 
Price 2003, Torii 2003) inside the SXC error circle. This is the brightest GRB ever detected 
by HETE-2 with the 30—400 keV fluence of 1.2 x 10~ 4 erg cm -2 , and precise and continuous 
follow-up observations were carried out by dozens of telescopes distributed around the world. 
Optical spectroscopic observations have determined its redshift as z = 0.1685 (Greiner et 
al. 2003), which is one of the closest ever known and is possibly related to the exceptional 
brightness of this afterglow. Moreover, spectra taken after several days reveal the evolution 
of broad peaks in the spectra characteristic of a supernova. The spectral similarity to SN 
1998bw (e.g., Galama et al. 1998, Iwamoto et al. 1998) and other energetic supernovae 
such as 1997ef provides strong evidence that GRB 030329 is associated with core-collapse 
supernovae (Dado, Dar & Rujula 2003, Hjorth et al. 2003, Stanek et al. 2003, Kawabata 
et al. 2003). In order to investigate the kinetic energy of GRB and the immediate vicinity 
of its progenitor, the early light curve is important. In this letter, we report the earliest 
detection of the optical afterglow of GRB 030329 starting 67 minutes after the burst. 

2. Observation and photometry 

Our observation was performed at Tokyo Institute of Technology using a 30cm-telescope 
(Meade LX-200) and an unfiltered CCD camera (Apogee AP6E) equipped with a front- 
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illuminated 1024x1024 CCD chip (Kodak KAF-1001E). The dark current images were sub- 
tracted from the obtained CCD images and then flat fielding was applied for all images. We 
used IRAF/noao/digiphot/apphot/phot packages to analyze the data. 

We started observing the preliminary SXC position at 12:44:13 UT on 29 March 2003, 
67 minutes (0.047 days) after the burst 4 . The magnitude of GRB afterglow at the very 
beginning was R c ~ 12.4 mag. This is the earliest detection of the afterglow of GRB 030329 
ever reported in literature. 

We continued observations for the rest of the night covering t ~ 0.05 — 0.30 days, and 
performed observations on the following two nights covering the period of t ~ 0.93 — 1.21 
days and t ~ 2.03 — 2.08 days, respectively, where t refers to the time since the burst onset. 
The exposure time of each CCD frame was 10 sec (0.047 < t < 1.21 days; when the after 
glow was relatively bright) or 30 sec (2.03 days< t). The magnitude of the comparison stars 
were calibrated using three stars in the same field of view, which has been calibrated in 
detail by Henden. Then we determined the magnitude of the OT in each frame relative to 
the weighted average of 15 bright comparison stars. 

Since the peak sensitivity of our camera is very close to the R c band, we have cali- 
brated magnitude of the OT by converting our magnitude system (Ri ns t) into the R c sys- 
tem assuming the color correction with the color (V — I) = 0.74 mag (Zeh et al. 2003) 
(a(r) = (Rinst) — 0.1514 x (V — I), where a(r) is the zero-point between the instrumental 
system and the standard magnitude). The statistical error is ~ 0.017 — 0.02 mag. Although 
the zero-point errors were found to be ~ 0.03 mag for the reference stars by Henden (2003), 
we found our data were consistently brighter than the i?-band filtered observations (e.g., 
Burenin et al. 2003). We therefore introduced additional zero-point correction (+0.11 mag) 
to match the light curves in the overlapped interval. The resulting light curve of the GRB 
030329 afterglow in the R c band is shown by combining data from other observations (see 
Fig 1). 

The light curve in the first day cannot be fitted with a single power-law function. We 
therefore tried to fit the light curve by two different forms of broken power-law functions. 
One is given by Beuermann et al.: F(t) oc [(t/t b ) a '^ n ' + {t/t b ) a ''z n ']~ l l n \ where t b is the break 
time, and n' provides a measure of the relative width and the smoothness of the break. The 
other is a "double-broken power-law" function with two breaks with the following form: 
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1 ) 1 i(t/t b2 r> n + (t/t b2 r^ 1/n (t>t hl ) 1 ) 

where t b i and t b 2 are the break times and n provides a measure of the relative width 
and the smoothness of the break. Here we excluded the "bump" at t ~ 0.08 — 0.09 days, 
which is discussed in section 3.3. 

We found that the former is not acceptable with a reduced x 2 of 1.72 (285 d.o.f) whereas 
the latter improves the fit significantly (the reduced x 2 1-06 with 283 d.o.f) (see Fig 2). 

As a result, it is well described by a broken power-law of the form; a\ = 0.88 ± 0.01 
(0.047 < t < t b i days), a 2 = 1.18 ± 0.01 (i bl < t < t h2 days), a 3 = 1.81 ± 0.04 (t h2 < t < 1.2 
days), where t b i ~ 0.26 days and t h2 ~ 0.54 days, respectively, and n = 18.8 ±5.1. Here, 
cki is determined by essentially the full Tokyo Tech data. The earliest phase of the light 
curve is well fitted by the single power-law with its index a\. a 2 and 0:3 are determined by 
measurements reported by Burenin et al. and GCN (see caption in Fig 1). 



3. Discussion 

3.1. Light curve at 0.05 < t < 0.26 days 

We have presented a light curve of the early phase of the optical afterglow of GRB 
030329 starting 67 minutes after the burst. This is the earliest detection of GRB 030329 
ever reported (Peterson & Price 2003, Torii et al. 2003, Uemura et al. 2003). 

Burenin et al. (2003) reported follow-up observations of GRB 030329 as early as 6 
hours after the burst, using BVRI filters. In each of the filters, they observed a gradual 
flux decay which can be accurately described as a power-law F v oc t~ L19 . They also found a 
characteristic break in the light curve t b rk ~ 0.57 day, after which the afterglow flux started 
to decline faster. The power-law slopes of the light curves changed from —1.19 to —1.9 for 
t > t br k- Notably, this break is nicely consistent with our second break (t b2 ) within error 
bars, where the power-law slopes changes from —1.18 ± 0.01 to —1.81 ± 0.04 after the break 
of t h2 ~ 0.54 days (see above). The first break (t bl ) is not discussed in Burenin et al, since 
they started their observations just around this break (t bi ~ 0.26 days). 

Price et al. (2003), Burenin et al. (2003), and Tiengo et al. (2003) found that the results 
of their observations are consistent with the model where the afterglow emission is generated 
during the deceleration of the ultra-relativistic collimated jet. They found that the break in 
the power-law light curve, at t ~ 0.5 — 0.6 days, can be interpreted as the jet-break, i.e., the 
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break which occurs when the angular structure of the ultra-relativistic collimated jet becomes 
observable (Sari, Piran & Halpern 1999, Rhoads 1999). This interpretation is supported by 
the fact that the break occurred simultaneously in different colors. Furthermore, a change in 
power-law slope from —1.19 to —1.9, is approximately consistent with that generally observed 
in jet-break. Therefore, our major concern in this paper is to understand the nature of the 
first break, tbi, and examine the consistency between the above scenario in the frame work 
of standard GRB fireball theories (Piran et al. 1999). 

3.2. Break at t ~ 0.26 days 

There are two possible break frequencies in the spectra, z/ m and z/ c , where z/ m is the 
synchrotron frequency, u c is the cooling frequency above which electrons lose their energy 
rapidly by synchrotron radiation (Sari, Piran & Narayan 1998). Since z/ m and u c are the 
functions of time, a break in the light curve could be observed when u c and/or z/ m crossed 
over the observed frequency z/r. Therefore, we examined six possible cases to understand 
the first break (tbi); according to the relation between z/ m , v c and z/r. In the standard 
GRB scenario, u c < z/ m is often called "fast cooling" since all electrons cool rapidly, whereas 
v c > v m is referred to "slow cooling" since only the high energy population of electrons 
cool efficiently We will also extend our discussion to discriminate between "a homogeneous 
interstellar medium (ISM) model" (e.g., Sari & Piran 1999) and "a pre-existing stellar wind 
model" (e.g., Chevalier & Li 1999) for the GRB environment. The relationships between 
observed spectral index and model predictions are compared in Table 1. 

We first consider the case where both u c and z/ m are above the observed optical frequen- 
cies (z/r < u m < u c or z/ R < u c < z/ m : case (3) and (6) in Table 1). In these two situations, 
observed flux at z/ R should increase with time, which strongly conflicts with the observed 
declining light curve. On the contrary, if the both cooling frequency v c and the minimum 
frequency u m are below the optical band (v c < z/ m < z/r or z/ m < z/ c < z/r: case (1) and (4) 
in Table 1), the predicted optical spectral index would be (3 — p/2, where p is the electron 
spectral index. Since the photon spectral index of this afterglow was /3 = 0.66 at t — 0.26 
days (Burenin et al. 2003), we expect p = 1.32, which is unusually flat for an electron 
population accelerated in a GRB. Furthermore, power-law index in the light curve should 
be a = ~ 0.49, which is again inconsistent with a± ~ 0.88 obtained with our data. 

Case (5) u c < z/r < z/ m in Table 1 is also ruled out since the predicted power-law 
index a = 0.25 (Sari, Piran & Narayan, 1998) is too flat compared to the observed value 
of ol\ ~ 0.88. Therefore, we argue that the possible solution is v m < u R < v c for the time 
region of t < tbi- In this case, however, if the burst occurred in pre-existing stellar wind, the 
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optical decay slope is predicted to be a = y - ~ 1.49, with 5 = 2 for a wind model 
(Panaitescu, Meszaros & Rees 1998), which is quite steeper than that observed, and hence 
we can rule out wind-interaction model. In summary, u m < z/r < u c and ISM model (case 
(2) in Table 1) is the only possible solution to reproduce both the temporal/spectral index 
of the optical afterglow of GRB 030329 at 0.05 < t < 0.26 days. 

In such a slow cooling case, time variation of afterglow flux is given by F oc t -3 ^ -1 )/ 4 
for v < v c and F oc t~( 3p ~ 2 )/ 4 for v c < is, respectively (Sari, Piran & Narayan 1998). By 
assuming ot\ = 0.88, the electron spectral index is estimated as p = 2.17. Note that, this 
electron spectral index agrees well with those of electrons accelerated in relativistic shock 
waves (e.g., Dado, Dar & Rujula 2003). Furthermore, we expect that power-law slope of 
the light curve would change from 0.88 to 1.13 for v c < v. Again, this is approximately 
consistent with the observed spectral index after t b i, where a 2 = 1.18 ± 0.01. Therefore, we 
conclude that the first break in the optical afterglow light curve at tbi is the most probably 
cooling-break where the cooling frequency crosses down the observed optical frequency. 

Under this assumption, we can determine important physical parameters for the GRB 
emission. For example, we can estimate e B and e e , where e# and e e are the fractions of 
the shock energy given to magnetic field and electrons at the shock (Sari, Piran & Narayan 
1998). In case of slow cooling, t m < £r < t c would be expected. Since we started our 
observation 0.047 days after the burst, we can limit the range of t m as t m < t i, s = 0.047 
days. For t c = t hl = 0.26 days, t m < 0.047 days, E = 10 52 erg, n = 1 cm" 3 , v = 0.5 x 10 15 
Hz, we obtain 

-<°H^)~W(!f)~W < 3 > 

We can also constrain the peak time of reverse shock days after the burst (Sari & 
Piran 1999). 



t A ~ 0.03 ( *L x ±) " ( ~ V^L x ^) " 2 (4) 
V0.05 0.1/ V 1.0/ V1.0 0.5 / y ! 

The Lorentz factor depends on time, 7(4) ~ (3 J B/2567rnm p c 5 t 3 ) 1/8 (Piran 1999). And 
the magnetic field is calculated using €3 and the Lorentz factor by B = •yc^32neBnm p . As- 
suming E = 10 52 erg, the Lorentz factor and the magnetic field strength at two characteristic 
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breaks time are 7 = 9.7, B = 0.86 gauss at t h i = 0.26 days and 7 = 7.4, B = 0.64 gauss at 
t h2 = 0.54 days, respectively. 

The values of e B and e e estimated in the preceding chapter are in good agreement with 
the averages of these parameters for GRBs calculated by Panaitescu & Kumar (2001), which 
are loge B = -2.4 ± 1.2 and e e = 0.062 ± 0.045. 

3.3. Bump at t ~ 0.08 - 0.09 days 

Finally, we comment on a small "bump" of the light curves at t ~ 0.08 — 0.09 days 
(Wnp) with an amplitude of ~ 0.1 mag. Uemura et al. (2003) reported a change of slopes 
from 0.74 to 0.95 days at t = 0.085 days. However our earliest data at t < 0.08 days has a 
slope steeper than 0.74. The light curve at t > 0.09 days lies on the extrapolation of this 
earliest segment. We consider this feature as a bump rather than a break. 

Short time variabilities, i.e., "bumps and wiggles", may be associated with the for- 
ward/reverse shock structures along the afterglow emitting regions (Kobayashi & Zhang 
2003), repeated energy injection from the central engine, or fluctuation in the density of the 
interstellar medium (Nakar, Piran & Granot 2003). 

First, we can rule out a case with the forward/reverse shock structure, since it predicts 
the light curve should not have the same power-law index before and after the bump. A case 
with repeated energy injection is also ruled out since after the injection the light curve after 
the bump should have the same power-law decay slope, but with a larger normalization. 
Therefore, we conclude that the bumps in the light curve is likely due to the fluctuation in 
the external density of the interstellar medium (Nakar, Piran & Granot 2003). 

We can estimate the distance from the central engine R(t) (Piran 1999) and the density 
variation (Nakar, Piran & Granot 2003) at tbump- 



(n/rao) ~ LI x [(F,/F )/l.l] 4 /( 1+p ) 
We find that density is enhanced about 10% at the distance of 2.2 x 10 17 cm. 



(6) 
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4. Summary 

We observed extremely bright optical afterglow of GRB 030329 from 67 minutes after 
the burst. Our observational results show that the shocked electrons are in the slow cooling 
regime with an electron index of 2.17 in this burst, and that the burst occurred in a uniform 
ISM, that is, GRB 030329 can be understood very well in the predicted "standard" model. 
We conclude the first break changes the power law index by ~ 0.3, consistent with the 
cooling-break in the frame work of the standard external shock model. 

We are grateful to the other members of the HETE-2 Ops Team for providing the 
location of GRB 030329. We are also grateful to R. Burenin and his collaborators for kindly 
providing us the numerical values of the filtered data and to K. Torii for giving useful advises 
on photometry. This work is supported by the Grants-in-aid for Scientific Research Program 
from the Ministry of Education, Science and Culture of Japan to NK. 
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Fig. 1. — Light curve of the optical afterglow of GRB 030329. The filled circles are our 
observations, the filled triangles come from Burenin et al. (2003), and the rest comes from 
GCN 2028, 2029, 2034, 2041, 2050, 2056, 2058, 2064, 2067, 2070, 2071, 2074, 2077, and 
KAIT. These magnitudes were translated using the standard sequence by Henden (2003). 



- 12 - 



■g 

'to 

CD 
DC 



CM 
=3 

■g 

'co 

CD 
DC 



U^J 

Reduced chi A 2 = 1 .72 (285 d.o.f) 




-», 



-I I L. 



Reduced chi A 2 = 1 .06 (283 d.o.f) 
' ■ i i i i ■ ■ ■ 




_i i i i i_ 



. 1 



Time (Days after the burst) 



Fig. 2. — Residuals of the light curve of the optical afterglow for the two models, top; the 
residual from a broken power-law which is containing a single, broad break, bottom; the 
residual from a double broken power-law (see text). 
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Table 1. Predicted decay slopes for various theoretical models 
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